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We study the far field spatial distribution of the quantum 
fluctuations in the transverse profile of the output light beam 
generated by a type II Optical Parametric Oscillator below 
threshold, including the effects of transverse walk-off. We 
study how quadrature field correlations depend on the polar- 
ization. We find spatial EPR entanglement in quadrature- 
polarization components: For the far field points not affected 
by walk-off there is almost complete noise suppression in the 
proper quadratures difference of any orthogonal polarization 
components. We show the entanglement of the state of sym- 
metric intense, or macroscopic, spatial light modes. We also 
investigate nonclassical polarization properties in terms of the 
Stokes operators. We find perfect correlations in all Stokes 
parameters measured in opposite far field points in the di- 
rection orthogonal to the walk-off, while locally the field is 
unpolarized and we find no polarization squeezing. 

PACS number(s): 42.50.Lc, 42.50.Dv, , 42.65. Sf, 42.50.Ct 



I. INTRODUCTION 

Polarization [1] and transverse spatial [2] degrees of 
freedom of light beams interacting with nonlinear media 
have been extensively studied in the last decade. The se- 
lection of special spatial modes of the transverse profile 
of a light beam down-converted by a quadratic crystal 
provides an interesting example of a polarization entan- 
gled state [3] . In this kind of experiments [3] the fluores- 
cence, or rate of photon pairs production, is low (single 
photon regime). Recently there has been an increasing 
interest for polarization entanglement in continuous vari- 
able regimes, where intense light beams, with high fluxes 
of photons, are generated. In this case the detection no 
longer resolves single photon events [4,5] . The interest in 
such macroscopic or multiphoton systems is partly due to 
possible applications of continuous variables in quantum 
communications [6], quantum information [7], mapping 
from light to atomic media [8], and quantum teleporta- 
tion [9] . Most works on continuous variable regimes [4-9] 
are concerned only with temporal features of light beams, 
while our aim in this paper is to study polarization en- 
tanglement between spatial modes of intense light beams, 
when intensities are continuous variables. 

Interesting polarization effects arise in type II phase 
matching, when a pump field is down-converted in a 



quadratic crystal in two orthogonally polarized fields [3] . 
Parametric down-conversion (PDC) can be increased us- 
ing an intense pump pulse or by means of a resonant 
optical cavity. The case of an intense pump is consid- 
ered in Ref. [10]. In this paper we study the situation of 
an optical cavity, that is a Optical Parametric Oscillator 
(OPO). In the OPO the fields resonate in a cavity and 
therefore an intense laser-like beam is down-converted 
above threshold. The cavity enhances the rate of produc- 
tion of photons pairs, increasing the gain and providing 
active filtering of the frequency bandwidth [11,12]. The 
spatial distribution of the quantum fluctuations close to 
threshold is dominated by weakly damped modes that be- 
come unstable at the threshold of the OPO, as discussed 
in Sect. II. This fluctuating spatial structure, known as 
'quantum image', has been extensively studied for type 
I OPO, where polarization does not play any important 
role [13-15]. Here we analyze the spatial quantum fluc- 
tuations for type II phase matching, considering the po- 
larization and the transverse spatial degrees of freedom. 
We also study the effects of the transverse walk-off be- 
tween the orthogonally polarized signal and idler down- 
converted fields. 

The characterization of the spatial and polarization 
properties of the down-converted light is given in two 
different ways, discussed in Sect. Ill and in Sect. IV, 
respectively. In Sect. Ill we study Einstein-Podolsky- 
Rosen (EPR) correlations [16] between polarization- 
quadratures components. A precedent of macroscopic 
EPR experiments in OPO are those of Ref. [17,18], but 
they do not refer to spatial EPR, distinguishing signal 
and idler by their polarization. Spatial EPR was theo- 
retically considered in Refs. [19,20] in type I-OPO. We 
build on these results by considering the polarization de- 
gree of freedom and including the walk-off effects in our 
treatment. Our main finding is that for the intersection 
points of the far field rings there is noise suppression 
much below the standard quantum limit in the proper 
quadratures combinations of any orthogonal polarization 
components of the critical modes. We show that the en- 
tanglement observed in single photon regime [3] survives 
for the multipohoton state generated by an OPO near 
below threshold. 

In Sect. IV we analyze the issue of macroscopic polar- 
ization entanglement in terms of Stokes operators. These 
operators are related to the intensity of different polar- 
ization components of light. There is a recent experi- 
mental observation of polarization squeezing in three of 
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the four Stokes parameters [21], while macroscopic po- 
larization entanglement in terms of Stokes parameters 
has been considered in Ref. [4] and demonstrated in Ref. 
[5]. The situations considered in Refs. [21,4,5] do not 
involve transverse spatial degrees of freedom. Quantum 
properties of Stokes operators taking into account the 
spatial transverse dependence of light seem to be first 
considered here and in Ref. [10]. We study both the 
local properties of Stokes parameters in the transverse 
far field profile of the down-converted beam (Sect. IV A) 
and the possibility of entanglement between spatial far 
field modes (Sect.IVB). As a main result we show that 
there exists macroscopic polarization entanglement be- 
tween the beams measured at the intersection points of 
the far field rings. We find perfect quantum correlation 
at zero-frequency in all Stokes operators. 



II. INPUT/OUTPUT RELATIONS AND FAR 
FIELD CHARACTERISTICS 

We consider a type II OPO below threshold, in the 
undepleted pump approximation. In this approximation 
the pump field is described by a fixed classical variable 
A . In the mean field approximation the signal (A\, or- 
dinary x polarized) and the idler {A 2 , extraordinary y 
polarized) fields obey the following Heisenberg operator 
equations [14]: 

<Mi = -71 (1 + *Ai - «aiV 2 + p 1 d y )A 1 + ^nA Q A\ 
+ y^A\ n (1) 
d t A 2 = -72(1 + iA 2 - ia 2 V 2 + p 2 d y )A 2 + ^yvnA A\ 

+ v/^Ar (2) 

where 7* are the cavity linewidths for the signal and idler 
fields, Aj are the cavity detunings, V 2 is the two di- 
mensional transverse Laplacian that models the effect of 
diffraction in the paraxial approximation, and at are the 
diffraction strengths. In type II phase matching a trans- 
verse walk-off arises between the signal and idler fields. 
It is described by the drift terms pid y Ai. Af 1 are the 
the input field operators, describing the fluctuations of 
the signal and idler modes entering through the partially 
transmitting cavity mirror. 
From the Fourier transform: 



where we have introduced the operator vectors 



d k f duj -,t s 



2tt J V2^ 
d 2 k f du 



i(k-x-U)t) 



A\x,t) = J ^ J ^e^-^AH-k, 



2tt 



(3) 
(4) 



we obtain the following algebraic linear relation, giving 
the intracavity mode operators C>(k, u>) in terms of the 
input fields v m (k, u)) 



0= ( Mk,w) 
yAt(-k,-u) 



and the matrices 

' 71(1 +iAi(fc,w)) 



r = 



with 



A?(k,w) 



AT\-k,-u;) 




72(1 - iA 2 (-k, -w)) 



Aj(k,u) = Aj + a 3 \k\ 2 + Pj ky - uj/jj. 



Using the input/output relations i> out 
we obtain the output fields: 



(6) 

(7) 

(8) 

(9) 
[22] 



A$(k, u) - u h2 (k, Lo)A\ n 2 (k, w) + v h2 (k, u)A 2 y(-k, -lo). (10) 



The coefficients of the input/output transformation are: 
2[l-iA 2 (-£,-w)] 



[1 + iAi(fc, w)][l - iA 2 {-k, -u)] - \A \ 2 

2A> 

[1 + iA 1 (k, w)][l - iA 2 (-k, -co)] - |A | 2 



-1 (11) 



(12) 



and U 2 , V 2 are obtained interchanging the indexes 1 and 
2 in Eqs. (11) and Eq. (12). 

Assuming that the input signal and idler fields are in 
the vacuum state, it is immediate to obtain the non- 
vanishing second order moments of the output fields 

(A™* \k,uj)A° 1 ut (k',(j')) = |U 1 (fc,w)| 2 <5(fc-fc')<5(^-'A (13) 

(A° ut \k,u)A° ut (k',Lu')) = \V 2 (k,uj)\ 2 8{k - k')5{u - w'), (14) 

(A° ut (k,uj)A 2 ut (k',Lu')) = C/i(fc,w)t/ 2 (-fc,-cj) (15) 
6(k + k')6(uj + uj 1 ) 

and the corresponding hermitian conjugate ones. From 
these moments, quadrature and intensity correlations can 
be analytically calculated for a transversally homoge- 
neous pump A n . These calculations can be simplified 
considering the unitary conditions (see, for instance, Ref. 
[14]) and the symmetries of the system. From Eqs. (13)- 
(14) we obtain the local Far Field (FF) intensity. From 
Eq.(15) it is evident that correlations between signal and 
idler are non vanishing only between symmetric points of 
the far field. Such two-mode correlations are due to the 
quadratic form of the Hamiltonian describing, within a 
linear approximation, the OPO below threshold. In the 
following, and in order to avoid unphysical divergences 
[23] , we consider that the fields are integrated over a de- 
tection region Rr of area a: 



LP = Tu l 



(5) 



A out {k) 



[ dk'i 



A out (k'). 



(16) 
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At the threshold for signal generation (A^' = 1), and 
for a negative total detuning (71A1+72A2 < 0), an insta- 
bility at finite wavenumber and with nonzero frequency 
appears, as extensively discussed in Refs. [24,25]. The 
modes that become unstable at threshold are determined 
by the relations Ai(fc, 0Jn(k)) = and A 2 (k, wjy(fc)) = 0, 
where 



uj H (k) 



7i72 
7i + 72 

[Ai - A 2 + (ai - a 2 ) 



(17) 



+ {Pi + P2)ky\ 



is the frequency that becomes undamped at threshold 
(Hopf bifurcation) . 

The unstable critical modes lie on two rings of the far 
field given by: 

71A1 +7 2 A 2 + (7iai +7 2 a 2 )|£| 2 ± (71P1 - l2P2)k y 

= 0. (18) 

If the relative walk-off 71P1 — -f 2 p 2 vanishes, Eq.(18) de- 
scribes a single far field ring. Therefore in absence of 
walk-off the signal and the idler far field distributions are 
superimposed, and an intense ring is observed. The two 
rings of Eq.(18) are clearly identified in Fig.l, where we 
represent the stationary mean intensity close to thresh- 
old: 



(A° ut "(k,t)Ar t (k,t) + A° 2 ut \k,t)A°2 Ut {k,t)) = 
£ f M\Vi(k,w)\ 2 + \V2(k,w)\ 2 ) 



(19) 



This figure is similar to the well-known experimental 
image obtained when there is no cavity, in spontaneous 
parametric down conversion [3,10]. We want to point 
out that the cavity introduces fundamental differences: 
in particular in the cavity case there is a threshold above 
which a pattern appears in the transverse profile of the 
fields. The modulus of the wave-vector of such pattern 
at threshold is identified in the noisy precursor below 
threshold [26]. The existence of a selected wave- number 
is an effect of the optical cavity, its value being deter- 
mined by the cavity detuning [24] . 

The main contribution to the integrals in Eq. (19) are 
the most intense frequencies components given by Eq. 
(17). Hence the intensity at the Hopf frequency 

{A° l ut \k,uj H {k))AT t (k : uj H (k)) + 
AT \k,w H (k))A° 2 ut (k,u; H (k))) 

is very similar to the one of Fig. 1. This picture of the FF 
shows that the intensity reaches the highest value at the 
intersection points of the two rings, where ordinary and 
extraordinary fields are superimposed. From Eq.(18) it 
is immediate to obtain the coordinates of these crossing 
points: 



±k H = (±K,0), with k% 



-71 Ai - 7 2 A 2 
7i<ii + 7 2 a 2 



(20) 



In our calculations we will also consider the FF modes 
on the rings for which the influence of the walk-off is 
stronger. These are the four points of intersection of the 
two rings and the y axis, with ordinates: 



±k c 1 

v± 2(7iai+7 2 a 2 ) 



[±(7ipi - 7 2/ o 2 ) 



(21) 



±\/ (liPi - J2P2) 2 - 4(7iAi + 7 2 A 2 )(7iai + 7 2 a 2 ) . 



We define the two external points by ±ky — (0, ±fc° + ), 
with ky + obtained from Eq. (21) with both + signs. 

For simplicity in the following we shall omit the label 
out, indicating with A\^ 2 the output fields, described by 
Eqs. (13)-(15). 



III. SPATIAL EPR ENTANGLEMENT BETWEEN 
QUADRATURE-POLARIZATION FIELD 
COMPONENTS 

The vectorial field is a superposition of linearly polar- 
ized components: 



A = Aie x + A 2 e y . 
By means of a wave retarder 



W = \ , ir 



1 

e l 



and a polarization rotator 

cos 6 sin 
— sin cos 



(22) 



(23) 



(24) 



we can obtain a field Are = WPA in any polarization 
state 

Are = (Ai cos + A 2 e tT sin Q)e x + 
(-Ax sin + A 2 e lT cos <d)e y . 

With a linear polarizer 

we can select a field polarization component, and, inte- 
grating over a detection region we obtain 

ire(M)= / dk'(A 1 (k , ,t)cosG + A 2 {k',t)e iT sm<d). 

By homodyne detection we can select a quadrature 
component of this polarization component. We define 
the quadrature ^ by 



A^it) = Are(ty* + (Are)Ht)e 



(26) 
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In any FF point, the arbitrary quadrature-polarization 
component (26) has vanishing mean value and a spectral 
variance which depends only on O, but it is independent 
of the choice of the phase factors and T. Integrating 
over a detection region R ± ^ of area a, much smaller of 
the variation scale of Ui and Vj, we obtain 

J dte™ (i^e (k, t)A* @ (k, 0)) = (27) 

<j[l + COS 2 QHV^k, LU\ 2 + \V!(k,-Lj\ 2 ) + 

sin 2 e(|y 2 (fc,^| 2 + |y 2 (fc, -uj)\% 

where a fixes the shot noise level. Therefore, in any 
far field point the variance of an arbitrary quadrature- 
polarization component is above the shot noise level a. 
We observe that if we place the detectors in the FF line 
k y = the variance results independent of all angles, 
including 6. We find 

J dte iut {Af e ((k x ,0),t)A^ @ ((k x ,0),0) = 

a[l + |^((fe,0),w| 2 + |Vi((fc x ,0),-u;)| 2 ], 

In other words, for k y — 0, the level of fluctuations is 
independent of the choice of the polarization state and 
quadrature, depending only on the position k x . 

Next, we consider the correlations between the field 
components detected from two spatially separated FF 
points. Such correlations arc not vanishing only for sym- 
metric FF points k and —k. We will show that the corre- 
lations between quadratures measured from these two po- 
sitions of the transverse plane show EPR entanglement. 
A two mode state is here defined to be EPR entangled if 
for two orthogonal quadratures Xi and Yi in each mode 
i (i = 1,2) the conditional variances V^mJj^-A^) and 
Kond(^il^2) are both less than 1, as discussed in Ref. [6]. 
The conditional variance is defined by 



V(A±gB) 
mm ; 

a V(Asn) 



(28) 



being V(A) the variance, and V(Asn) the shot noise 
level. The factor g is introduced to optimize noise re- 
duction and is experimentally obtained by an attenuator 
and a delay line [20]. We note that 



V(Xi - X 2 ) < V(X 1>S n) and 
V(Yi + Y 2 ) < V(Y hSN ) 



(29) 



is a sufficient condition for EPR entanglement, corre- 
sponding to the choice g = 1. The definition of EPR 
entanglement used here [6] provides a sufficient condi- 
tion for the inseparability criterion recently discussed for 
continuous variables systems in Ref. [28]. 

For a single-mode type II OPO below threshold, in 
which transverse effects are not considered, EPR corre- 
lations between signal and idler modes of different polar- 
izations have been predicted and experimentally demon- 
strated [27] . Recent investigations show the possibility of 



EPR entanglement between spatial regions of the trans- 
verse profile of the signal field of a degenerate optical 
parametric oscillator [19,20]. For type II phase match- 
ing we can consider two symmetric far field modes with 
x and y polarizations, respectively. Neglecting walk-off 
effects we would then find results equivalent to the degen- 
erate case in type I phase matching [19,20]. In addition 
of considering the effect of the walk-off, the novelty of our 
discussion here for type II OPO is that we can also con- 
sider how these correlations change with the polarization 
state. 

As an indicator to look for EPR entanglement in our 
case, we introduce the spectral variance of the difference 
of quadratures 

V g (±k,uj;$)= J dte iut {(Af @ (k,t)-g*A^ @ ,(-k,t)) (30) 

(i* e (fc,0)-. g i*/ e ,(-fc,0))), 

where the vector $ = (r, 9, T', 6', is the set of pa- 
rameters determining the polarizations (lower labels) and 
quadratures (upper labels) of the symmetric FF modes 
under consideration. The value of g giving the minimum 
variance (30) is 



Jdte^(A^ e ,(-k,t)Af e (k,0)) 
9 J dte^{Af l& , {-%, t)Af l@l {-k, 0)> 

From the output moments Eqs. (13-15) we obtain 



(31) 



V 3 (±fc,w;$) =cr (32) 
■| e i(*+*'+r') cos8[/l (fc ; w ) _ 5 * s in0'T/ 2 *(-fc, -co)| 2 

+ |e l (* + *'+ r '). 9 *sinO'[/ 1 (fc,-co)-coseK 2 *(-fc,w)| 2 
+ | e *(*+*'+ r ) sinei/i(-ft,-w)-s* cose'y 2 *(fc,w)| 2 
+ | e J (*+*'+ r ) g* cos O'C/i (- jfe, u) - sin OV 2 * (k, -lo)\ 2 



with shot noise level <r(l + |.g| 2 ). A variance below this 
shot noise level is a signature of squeezing. We are 
looking for the more stringent conditions (28) discussed 
above: EPR entanglement imposes the requirement that 
V g (±k, ui; <E>) goes below the shot noise level of Af e (that 
is we should find V g (±k, w;<&) < a), simultaneously for 
two combinations of orthogonal quadratures [6] . 

The general result (32) depends on many parameters. 
However it is important to note that V g only depends on 
the phases ^,^',r, T' through the independent combi- 
nations (* + *' + T) and (* + *' + T'). The dependence 
on the sum of quadratures angles ^ and is well known 
in other contexts. It is easily understood taking into ac- 
count that measuring by a single homodyne detector the 
noise in a quadrature of the difference of the spatial 

±k and polarization (r, 6, T', 6') modes is equivalent to 
the noise measurement described by Eq. (30) [27] . The re- 
sult that V g depends on independent variations of T and 
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T' only through their sum with the sum of the quadra- 
ture phases \& + means that it is equivalent to vary the 
selected quadrature changing the phase of the local oscil- 
lator or to shift both signal and idler fields by a proper 
phase with the phase retarders. 

In order to study how correlations change in differ- 
ent spatial points of the FF we consider two detection 
schemes, with detectors in symmetric FF points, which 
represent possible extreme cases: 

• In the first detection scheme the field is detected in 
the crossing points ±kn of the signal and idler rings 
in the FF, as represented in Fig. 2. Being these 
points in the line k y = 0, they are not affected by 
the transverse walk-off. In Fig. 2 we indicate the 
case in which in one point the polarization compo- 
nent x is selected and in the symmetric point the 
polarization y is selected. In Sect. Ill A we will 
show that for these special points of the FF any 
change in the selection of the polarization does not 
influence the results given by Eq. (32). 

• In the second arrangement, shown in Fig. 3, the de- 
tectors are located in a couple of symmetric points 
±fcy ° n the line k x = 0, where the walk-off effect 
is more pronounced. In Sect. IIIB we will con- 
sider the effect of changing the polarization state 
selected. We can distinguish two extreme cases: 
Fig. 3a shows the arrangement in which the most 
intense polarization components on the rings are 
detected. We name this scheme as "vertical bright 
scheme" . Fig. 3b shows the case in which the y (x) 
polarization component is selected on the intense 
lower (upper) x (y) polarized ring. We name this 
scheme as "vertical dark scheme" . 



A. EPR between far field modes unaffected by 
walk-off 



In the line k y = there is no effect of the transverse 
walk-off and the coefficients given by Eqs.(9), (11) and 
(12) have the following reflection symmetry: 



Uj{h X: ky — 0) 

Vj(k X ,ky=0) 



Vj(- 



= 0) 
0) 



(33) 



with j — 1,2. The results presented in this section are 
strongly dependent on this symmetry, generally present 
also in previous treatments of spatial EPR correlations 
and squeezing. 

Even with this symmetric form of the coefficients the 
variance given by Eq.(32) is a complicated function of 
many parameters. For the sake of simplicity we first con- 
sider Eq.(32) in the case of g = 1 (see Eq.(29)). The 
microscopic process of generation of twin photons with 
linear orthogonal polarization (6 = and 6' = 7r/2) 
suggests a natural choice for the relative phases of the 



polarizers. Hence we consider a case in which the polar- 
izers in the symmetric points ±kn have a relative phase 
fixed by 



e' = e + 7r/2. 



(34) 



Using this relation between 6' and 6, Eq.(32) becomes 
independent of 9 for g = 1 when the additional choice 



r' = r + 7T 



(35) 



is also made. With these particular choice of parameters, 
Eq. (32) for the ±kjj points (20) reduces to: 

v g=1 (±k H , w, [r, e, *, r + tt, e + tt/2, *']) = 

\e i (*+*'+Vu 1 (k H ,cj) + V 2 *(k H ,-uj)\ 2 
+ \e^+*'+ T )U l {k H ,-u J ) + V*{k H ,u J )\ 2 \ . (36) 

Eq.(36) explicitly shows that the fluctuations in any 
quadrature of the difference of symmetric spatial modes 
±kn, with relative polarizations fixed by Eqs.(34)-(35), 
are independent of the choice of the polarization reference 
9. 

With the above motivation for the relations between 
phase parameters, we analyze the EPR correlations in the 
±ku points (Fig. 2) fixing the parameters 9 = 0, T = n, 
9' = tt/2, T' = 2tt, so that i* e = Af and i*,' e , = if, 
and varying the quadrature angles "J/, \E ,/ . There is no 
loss of generality in this choice of T and 9 since Eq.(36) 
is independent of 9, and the phase T can be absorbed 
in the quadrature angles IP + As discussed previ- 
ously (see Eq.(29)), EPR entanglement is guaranteed - 
for some quadratures (<J/ + - if both the 'position' 
and 'momentum' operators 

A*{k H )-Af{-k H ) = (37) 
[A^kny* + H.c] - [A 2 {-k H )e^' + H.c] 
A^' 2 {k H ) + At + *l\-k H )= (38) 
[iiitfcffje** + H.c] + [iA 2 (-k H )e™' + H.c] 

show simultaneously a variance below the reference value 
a. This value corresponds to the shot noise level of 
both Af(k H ) and A* +7T/2 (k H ). Using Eq. (32) it turns 
out that (37) and (38) have the same spectral variance. 
Therefore, in the following we only consider Eq. (32) 
normalized to the shot noise level a for the 'position' 
operator Eq. (37), which is identified by the angles 
I = (tt, 0,tf, 0,7r/2, *')■ 

In Fig. 4 we represent the normalized variance 
Vg = i(±fcff, u>; for the points ±fc#. When this nor- 
malized variance is less than 1, we find EPR entangle- 
ment. We see a maximum noise reduction for fy + ty' = 0, 
and for uj = 0. Fig. 5 shows a cut of Fig. 4 for ^ + ty' = 
(continuous black line). A variation in the quadrature 
angle results in a mixing of squeezed and unsqueezed 
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quadratures degrading the entanglement: strongest de- 
grading effects are evident at small frequency (dashed 
black in Fig. 5). 

The fact that there is maximum noise reduction for 

+ \E r ' = can be understood by analogy with the re- 
sult for the single-mode non-degenerate OPO discussed 
in Ref. [27]. In fact, for the case of vanishing detunings 
of the signal and idler and for real pump, which was the 
one considered in Ref. [27], the squeezing direction cor- 
responds to the real quadratures = *S?' = 0). In our 
case we find an analogous result because we are also con- 
sidering a real pump, and the effective detunings, given 
by (9), also vanish: Aj(k — ±fc#,w = 0) = 0. 

Our calculations allows us to search for EPR entan- 
glement considering any couple of symmetric FF modes 
(distinguished by their transversal wavevector) and se- 
lecting any polarization. We have found EPR entangle- 
ment for the modes ikn for any choice of the polariza- 
tion component in — kji (varying 9), if in +kn we select 
the orthogonal polarization (9' = O + tt/2). In these far 
field points ±Jch, which are not affected by the walk-off, 
any mixing of the the signal and idler fields detected in a 
point is entangled with the field in the symmetric point, 
if this is also properly mixed (Eqs.(34)-(35)). 

So far we have considered the case g = 1 and we have 
found that a sufficient condition to guarantee EPR entan- 
glement between the modes ±kn is fulfilled (see Eq.29). 
If we now optimize the noise reduction considering g = g 
we obtain the results shown in Fig. 6. Comparing the 
normalized variance V g= \ (Fig. 4) with Vg (Fig. 6) we 
observe that with the choice g = g, EPR entanglement 
is observed in a larger frequency bandwidth. For a more 
direct comparison of the variances obtained for g = g 
and g =1, both quantities are plotted for \I> + = in 
Fig. 5: for small frequencies the results are very similar, 
while for increasing values of the frequencies the choice 
g — g allows to observe EPR entanglement even when 
this effect is lost for g = 1. 

Our strong EPR correlations have been obtained for 
O' = 6 + 7r/2. As already mentioned, this phase relation 
corresponds to the underlying process governing the cre- 
ation of twin photons with orthogonal polarizations. In 
addition, for j = lwe have obtained that the EPR corre- 
lations are independent of 6. If we decrease the relative 
angle between O and 9', we observe a progressive reduc- 
tion of the entanglement. In the limiting case 9 = 9' 
the correlation between symmetric spatial modes is van- 
ishing and V g > cr. In fact in this case no twin photons 
are detected. 



B. EPR between far field modes in the walk-off 
direction 

In this section we study possible EPR entanglement 
for symmetric points of the far field along the walk-off 
direction (y), in the arrangements shown in Fig. 3. An 



important effect of the transverse walk-off is breaking 
the reflection symmetry in the far field. This symmetry 
is generally broken for k y ^ 0. For the points shown in 
Fig.3 

Uj(k x = 0, ky) ^ Uj(k X = 0, -ky) J = 1, 2 (39) 

Vj(k x = 0, k y ) ± V 3 (k x = 0, -ky) j = 1, 2. (40) 

Following the considerations in the previous Section, 
we will also consider here the case of phase polarizations 
9, 9', T, T' fixed by Eqs.(34)-(35). For this special choice 
of phase relations, the variance (32) for arbitrary g re- 
duces to: 

V g (±k,Lu;$) = (41) 
(t{cos 2 9[|e j( *+*' +r) C/i(fc, u) + g*V 2 *{-k, -lo)\ 2 
+ \e^' + ^g*U 1 (k,-u J ) + V 2 *(-k,u;)\ 2 } 
+ sin 2 e[|e < <*+*' +r > L7i (-£, -w) + g*V* (jfc, u) | 2 
+ \e^'+^g*U 1 (-k,cu) + V 2 *(k, -c)| 2 ]}. 

The lack of reflection symmetry implies that the variance 
Vg{±k,u;§) depends now on the angle 9 and a sim- 
ple result analogous to Eq. (36) cannot be obtained when 
k y ^ 0. Eq.(41) depends on the sum * + + T: there- 
fore, without loss of generality, we can absorb the effect 
of the wave retarder in the phase of the local oscillator, 
fixing the angle r = tt as in Sect. Ill A. In the following 
we consider the dependence of Eq.(41) on the angle 9. 
Varying 9 different polarization components are selected 
locally. We will see that the selection of different values 
of 9 can improve or degrade EPR entanglement. In par- 
ticular we consider two values for the angle 9 (= 0, tt/2) 
leading to very different situations. 

First we consider the case represented in Fig. 3a, in 
which the polarizer is oriented so that the most intense 
linear polarization component is selected locally. We are 
selecting two critical spatial modes A\(— ky) and A 2 (ky), 
whose quantum fluctuations are weakly damped. This is 
the "vertical bright" detection scheme. In order to detect 
the intense polarization components at ±fcy the phase 9 
must be fixed at 

9 = tt/2 (42) 

and, given Eq.(34), 9' = it. Therefore the phases to be 
considered are $ = (tt, tt/2, ^, 0, tt, $'). We first consider 
the case g = 1. We look for EPR entanglement between 
the 'position' and 'momentum' quadratures 

A?(-kv)-A*(kv), if' +w/2 (-M + i* +T/2 (M- (43) 

We will only show the variance of the 'position' quadra- 
ture, the orthogonal one being equivalent. We obtain 
EPR entanglement, as shown in Fig. 7. Maximum noise 
reduction is obtained for ^ + = 0. The variance V g= \ 
normalized to a and for ^ + = is represented as a 



G 



function of the frequency lo in Fig. 8: the best entangle- 
ment is observed for lo = 0. Also in this case a varia- 
tion in the quadrature angle '3/ + results in a mixing 
of squeezed and unsqueezed quadratures degrading the 
entanglement (Fig. 8). We note an important difference 
with respect to the case for ±kn- in Fig. 5 the largest 
degradation of entanglement for ^ + ty' ^ was ob- 
served for vanishing frequency {lo ~ 0), while in Fig. 8 
we see that for vanishing frequency the entanglement is 
only partially degraded. The largest degradation occurs 
now for lo ~ ±0.4045 (two peaks in Fig. 8). This value of 
frequency coincides with the Hopf frequency loh for the 
mode ky, as can be easily checked by Eq.(17). 

In Fig. 9 we show the optimized (g = g) variance V g=g 
(41) normalized to the shot noise a. We obtain EPR 
maximum entanglement for + vp' = and for small fre- 
quencies. Both the minimum and maximum fluctuations 
are obtained in a bandwidth of frequencies centered in 
zero, as in the case for the points ±fcjy. The effects of 
the Hopf frequency found for g = 1 disappear for g = g. 
A clear difference with respect to the variance V g=g ob- 
tained in Sect. Ill A is the level of maximum noise sup- 
pression reached. Even if in both cases we observe strong 
EPR entanglement, in this "vertical bright" arrangement 
the quadratures correlations are reduced, as can be seen 
comparing Fig. 5 and Fig. 8. This reduction is caused by 
the walk-off. The fields in the critical modes ±ku -not af- 
fected by walk-off- have vanishing effective detunings (9) 
for the threshold frequency lor^h) = 0; while the fields 
in the critical modes ±fcy -in the walk-off direction- have 
vanishing effective detunings (9) for the threshold fre- 
quency Wff(±fcy) 7^ 0. The fact that the detunings do 
not vanish for lo = seems to be the mechanism re- 
sponsible for the reduction of squeezing at u> = for 
the modes ±ky On the other hand, also the unsqueezed 
quadrature is influenced by this effect, showing a reduced 
amplification with respect to the values obtained for the 
points ±fcjj. Comparing Figs. 6 and 9, we see that in this 
"vertical bright" arrangement the variance is less sensi- 
ble to deviations of the quadrature phases selected, with 
respect to the optimum choice ^ + = 0. In Fig. 9 we 
observe a broad interval of phases ^ + giving EPR 
entanglement for lo = 0. 

Next we consider the detection scheme of Fig. 3b. In 
this case the phase of the polarizer at +ky is fixed at 

6 = 0, (44) 

and $ = (7r, 0,^,0, 7r/2,*'). With this selection of the 
phases of the polarizers the intense field component arc 
filtered out. In this "vertical dark" detection scheme, the 
detected modes {A±{kv) and A 2 {— ky)) have low intensi- 
ties. The main point is that now we are considering non 
critical modes that are strongly damped at any frequency 
(see Sect. II). We evaluate again the spectral variances of 
the position and momentum quadratures. In this case 
the results obtained for g = 1 and g = g are completely 



different. We start considering Eq.(41) for g = 1. We ob- 
tain that Vg—i is always larger than a, therefore no EPR 
entanglement is observed (see Fig. 10). In the same way 
that in the vertical bright scheme, the largest fluctua- 
tions are observed at lo ^ 0. Fig. 11 shows a cut of Fig. 10 
for the quadrature for which there is a minimum in 
the direction * + 

We now consider the variance (41) for the best choice 
g = g. We find that V^, represented in Fig. 12, is re- 
duced below the shot noise level a for a large region of 
parameters. Therefore, with a proper choice of g, EPR 
entanglement is obtained also in this case. From Fig. 12 
we also see that for ^ + ^E 1 ' = only small entanglement 
would be observed, and in the region of large frequencies. 
In fact, the quadrature at which strong EPR effects are 
observed is * + = § ^ 0. Fig. 11 shows V g for this 
choice of ^ + Changing the walk-off parameter we 
have found that \& increases with the walk-off. 

In Fig. 13 we show a comparison of the best EPR en- 
tanglement (g = g, optimum ^ + ty') found for the three 
detection schemes considered, namely Fig. 2, Fig. 3a and 
Fig. 3b. We observe that the correlations are less impor- 
tant when we move from the detection scheme of Fig. 2 to 
the the "vertical bright" and "vertical dark" schemes of 
Fig. 3a and Fig. 3b. We also show the effect of increasing 
the walk-off: we can see that EPR entanglement in the 
"vertical dark" and "bright" schemes get worst increas- 
ing the walk-off strength. Obviously the results for the 
detection scheme of Fig. 2 are not influenced by walk-off. 

IV. STOKES OPERATORS 

In the previous Section we have seen how the selec- 
tion of different polarization components in the far field 
influences the quadratures EPR entanglement between 
symmetric FF points. The results we have obtained also 
show the effects of the transverse walk-off. In this Sec- 
tion our aim is to characterize the polarization properties 
of a type II OPO, when transverse walk-off is taken into 
account. The polarization state of the field in any point 
of the transverse plane can be characterized in terms of 
the Stokes parameters [33]. An operational definition 
of the Stokes parameters can be given using polarizers 
and retarders [33,34]. In the quantum formalism there 
are different ways to describe the polarization giving the 
same classical limit [35,36]. Here we consider the quan- 
tum Stokes operators Sj(k, t) (j = 0, 1, 2, 3), for each FF 
mode k, obtained replacing by creation and annihilation 
operators the corresponding observables in the classical 
definitions (see Ref. [37] and references therein): 

S {k,t) = iJ(M)ii(M) + i£(M)i 2 (M) (45) 
is the total intensity operator, 

&(M) = it(M)ii(M) - Al(k,t)A 2 (k,t) (46) 
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gives the difference between the x and y linear polariza- 
tions, 



S 2 (k,t) = A\(k,t)A 2 (k,t) + A\(k,t)Ai(k,t) (47) 

gives the difference between the 45° and 135° linear po- 
larizations, and 

S 3 (k,t) = -i(A\(k,t)A 2 (k,t) - A\(k,t)Ai(k,t)) (48) 

gives the difference between the right-handed and the left 
handed circular polarizations components [42] . The defi- 
nitions (45)- (48) correspond, except for a constant, to the 
the Schwinger transformation of the modes Ai(k,t) and 
A 2 (k,t) giving operators satisfying angular momentum 
commutation relations [38]: 



S (k,t),Sj(k',t')\ =0 (49) 
Sj(k,t),S k (k',t')] = 2ie jkl S l (k,t)8(t-t')a5 m , (50) 



with j, k, I = 1,2,3 [43]. The precision of simultaneous 
measurements of the Stokes operators is limited by the 
Heisenberg principle. For instance we have 



A 2 Si(M)A 2 S 2 (M') > -\([Si(k,t),S 2 (k,t')})\ 2 
= \(S 3 (k,t))\ 2 6(t-t'). 



(51) 



The Stokes vector S = (Si, S 2 , S3) can be represented 
in a quantum Poincare sphere, with a radius defined by 
(So) (see for instance Ref. [4]). Given the fluctuations 
of Si, the quantum states are not defined by points on 
the surface of this sphere, but rather they are defined 
by different volumes, such as spheres (coherent states) or 
ellipsoids (squeezed states). These quantum uncertainty 
volumes on the Poincare sphere have been confirmed by 
recent experiments [21]. The transformations Eqs.(23) 
and (24) introduced in the previous Section can be also 
visualized in the Poincare sphere. In fact they corre- 
spond to a rotation in the Poincare Sphere of an angle 
20 around the S3 axis and of — r around the Si axis. 



A. Far field local properties 

Given Eqs. (13-15) we obtain the stationary value of 
the average of the Stokes parameters. The average of So 
is given in Eq. (19) and for the others Stokes parameters 
we find: 

(Si(M)> = ^J M\Vi{k,u)\ 2 - \V 2 (k,Lu)\ 2 ) (52) 

(S 2 (k,t)) = (S 3 (k,t))=0. (53) 

The average of S 2 and S3 vanishes in any point of the 
far field, because any signal or idler photon has the same 



probability to be measured along the 45° and 135° po- 
larizations directions. The same is true for the left and 
right circular polarizations. The equal average intensi- 
ties at the output of the beam splitters are subtracted, 
giving vanishing values of (S 2 ) and (S3) independently of 
the relative intensity of the signal and the idler. In Fig. 
14 we show the far field spatial profile of (Si): the lower 
ring is dominated by the linear polarization x while the 
upper one is dominated by the y polarization. If there 
was no walk-off, (Si) would vanish in all the FF, while 
in our case (pi,p 2 7^ 0) it only vanishes along the di- 
rection k y = 0. Therefore, in the ±fc# points we have 
an intense field (see So in Fig. 1) with vanishing aver- 
age of the Stokes vector S and with variances of Sj not 
limited by the Heisenberg principle, since from Eq. (51) 
\([S i (k,t),S j (k',t')})\=0,i^j. 

The parameter that corresponds to the classical 
characterization of the polarization state of a quasi- 
monochromatic field, is the second order polarization de- 
gree 



£j=ito(M)) 2 
(S (k,t)) 



(54) 



varying from P 2 = for unpolarized light, to P 2 = 1 for 
completely polarized light. In Fig. 15 we observe how 
the polarization degree, that reduces to 



Pi = 



\{Si)\ 
(So) 



(55) 



varies in the FF: In particular, the intense FF rings are 
always polarized except around the line k y = 0, where P 2 
vanishes. Therefore, the field in the points ±fcy is un- 
polarized in the ordinary sense. However the concept of 
polarized and unpolarized light needs to be generalized 
in quantum optics [30,37-40]. The fact that (Sj) = 
(j = 1,2,3), so that P 2 vanishes, does not guarantee to 
have an unpolarized state from a quantum point of view. 
Rather one has to consider the values of the higher or- 
der input moments of Sj. In Ref. [30] it was shown that 
in the single-mode type II PDC, the squeezed vacuum is 
unpolarized in the ordinary sense (P 2 = 0) due to the dif- 
fusion in the difference of the signal and idler phases. On 
the other hand, due to the twin photon creation, there is 
complete noise suppression in the intensity difference of 
the two linearly polarized modes, i.e. in Si, leading to 
polarization squeezing [29]. Due to this anisotropic dis- 
tribution of the fluctuations in the Stokes vector S there 
is a "hidden polarization" [30], which has been observed 
experimentally recently [37] . We note that similar hidden 
polarization should be observed locally in the transverse 
near field plane of a type II OPO. Our interest here is 
in the far field plane, where twin photons are spatially 
separated. 

In order to characterize the far field polarization prop- 
erties, we proceed to evaluate the variance of the Stokes 
parameters. We define the spectral correlation function 
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Ti(k,k',n)= (56) 
j dt^Siik^Siik 1 ,0)) - (^(M))(^(fc',0)}]. 

Given the moments (13-15) and using the moments theo- 
rem [41], we obtain non vanishing contributions only for 
k! = k (self correlation) and for k' = —k (twin photons 
correlation). For O = 0, corresponding to integrating 
over a time interval long enough with respect to the cav- 
ity lifetime, the self correlation of Si is 

ri(M,o)= (57) 

v 2 J ^[|C/i(fc,o;)| 2 |yi(fc,a;)| 2 + |C/2(fc,w)| 2 |V 2 (fc, W )| 2 ] 

while the twin photons correlation for symmetric FF 
points is 

T 1 (k,-k,0) = -T 1 (k,k,0). (58) 

It can be easily shown that the second and the third 
Stokes operators have equivalent variances: 

T 2 (k,k',uj) =T 3 {k,k',uj). (59) 

For the self correlations we obtain 

r 2 , 3 (M,o)= (60) 
v 2 J ^[|C/ 2 (fc,o;)| 2 |^(fc, W )| 2 + |C/i(fc,a;)| 2 |y2(fc,w)| 2 ] 

and for the twin correlations 

r 2 ,3(A,-M)= (6i) 

ff2 / Tj^Pi&uM-k' -u)V 2 (k,w) + c.c] 

For symmetry reasons Eq.(57) and Eq.(60) give iden- 
tical results in the crossing points of the rings k = kn 
(see Fig. 16). Since the fluctuations are isotropically 
distributed in S no "hidden polarization" is observed in 
these points. Let us now consider the possibility of quan- 
tum effects. First, the shot noise level for all the Stokes 
parameters (i.e. k, 0) evaluated on coherent states) 

is given by the average total intensity (So(k, t))a (see Eq. 
(19)) [4,30]. Inspection of Eq. (57) and Eq. (60) shows that 
all the Stokes operators have classical statistics in any FF 
point, as shown in Fig. 16. In these diagrams we rep- 
resent the normal ordered variances. These are obtained 
from the variances Eq.(57) and Eq.(60) after subtraction 
of the corresponding shot noise. We see that these are 
positive quantities in all the far field, i.e. no polarization 
squeezing appears [29] . The physical reason is that twin 
photons are emitted with symmetric wave- vectors (sym- 
metric FF points), while locally no correlations between 
orthogonal polarizations are observed. This motivates us 
to consider in the next Section the correlations between 
the Stokes operators of twin beams. 



In conclusion, when there is walk-off the polarization 
state (P 2 ) varies in the FF, the fluctuations in the Stokes 
parameters are above the classical level in all the far field 
and no polarization squeezing is observed. For k y = 

and in particular for ±fc#, P 2 = and fluctuations are 
isotropically distributed in all Stokes parameters so that 
the field is completely unpolarized. This result would 
apply to all the FF plane if there was no walk-off. 

B. Far field correlations 

From Eq. (58) and Eq. (61) we see that Si for sym- 
metric FF modes is anticorrelated, while 5 2 and S3 are 
positively correlated. The physical reason for the sign of 
these correlations is always the underlying twin photon 
process which creates pairs of photons with symmetric 
wave-vector and orthogonal polarizations x and y, lead- 
ing to a positive correlation of the corresponding beam 
intensities. 

These considerations suggest to look for noise suppres- 
sion in the following superpositions of Stokes operators 

Di{k,k',w) = J dte iajt {([Si{k,t) + Si{-k',t)] (62) 
[Si(k,0) + Si(-k',0)})} 

and 

Di(k,k',uj) = J dte iut {([Si{k,t)- Si(k',t)] (63) 

[&(M) -&(#,())]>} 

with i = 0,2,3. We note that Di(k,k' ,u) — D (k,k' ,uj) 
(see the definitions Eq. (45) and (46)). In addition, 
it follows from Eqs.(60)-(61) that D 2 = D 3 . These 
quantities are non-vanishing for symmetric points k' — 
—k. From the unitarity of the transformation (10) we 
know that \Ui{k,uj)\ 2 = \U 2 (-k, -cu)\ 2 and \Vi(k,u)\ 2 = 
\V2(— k, — oj)\ 2 , so that from Eqs. (52) and (53) we obtain 
that 

(Si(k,t) + Si(k',t)) = (S i (k,t)-S i (k',t))=0 (64) 

(i = 0,2,3). Therefore Eqs. (62) and (63) actually define 
variances. 

From Eq. (57-58) we obtain 

Di(k, -k, 0) = Ti(k, k, 0) + Ti(-k, -k, 0) + (65) 

Ti(-k, £,0) + ri(fc,-fc,o) = 

for any k: There is complete noise suppression in the 
sum of the Stokes parameters Si evaluated in symmetric 
regions of the transversal field. Therefore the normal 
ordered variance is equal to minus the shot noise taking 
nonclassical negative values. In conclusion, due to the 
twin beams intensity correlations, we find entanglement 
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for Si evaluated in any symmetric FF points and for any 
pump intensity. 

The Stokes operators generally cannot be simultane- 
ously measured with infinite precision (sec Eq. (51)). 
Also the superpositions of Stokes operators involved in 
(62) and (63) are in principle limited by the Heisenberg 
relations. However we have that 

([S 1 (fc,i) + 5i(-fc,t),S 2 (fc,0-^2(-fc,i')]> - (66) 
2ia(S 3 (k, t) - S 3 (-k, t))6(t -t')=0 

as follows from Eq. (64). Therefore there are super- 
positions of the Stokes operators whose measurement is 
not limited by Heisenberg relations because the average 
of their commutator vanishes. In the same way, also the 
other superpositions of Stokes operators for any couple of 
symmetric points can be simultaneously measured with 
total precision. This result opens the possibility to ob- 
serve noise suppression not only in the first Stokes oper- 
ator superposition (65), but also in the other superposi- 
tions Eqs. (63). However, we obtain that 

D 2 (k, -k, 0) = r 2 (fc, k, 0) + T 2 (-fc, -k, 0) - (67) 

r 2 (-fc, fc,o)-r 2 (fc,-M)^o, 

so that D 2 (k, — k, 0) generally does not vanish for arbi- 
trary k. The profile of the normal ordered D 2 (k, —k, 0) 
is represented in Fig. 17: we observe that in most part 
of the far field where there is large light intensity this 
quantity is positive, giving classical statistics. However 
there is a bandwidth of small wave-vectors k y ~ for 
which the normal ordered D 2 is negative and therefore 
quantum effects are observed. This small bandwidth be- 
comes smaller when increasing the walk-off, as can be 
seen comparing Fig. 17 -obtained with p 2 = 1- with 
Fig. 18 -obtained with p 2 = 0.5-. 

In particular, along the k y = line, we obtain that 
D 2 (k x , — k x , 0) = 0. This result easily follows from the 
symmetry Eq. (33). Therefore, along the k y = line 
we have D\ = D 2 = D 3 = which indicates perfect 
polarization entanglement between symmetric FF modes. 
In summary, in the direction orthogonal to the walk-off 
(k y — 0), we show complete noise suppression in properly 
chosen symmetric modes superpositions of all the Stokes 
parameters. Along this line the two points ±fc# are of 
special interest because they have a large photon number. 

We finally point out that, the situation studied here is 
different from the case of bright squeezed light considered 
in Ref. [4] and that the perfect correlations obtained be- 
tween Stokes parameters cannot be used here to obtain 
an EPR paradox: in fact from Eq. (51) and Eqs. (52- 
53) we obtain that the Heisenberg principle imposes no 
limits in the local variances of the Stokes parameters. 



V. CONCLUSIONS 

We have investigated the EPR entanglement between 
quadrature-polarization components of the signal and 
idler fields in symmetric FF points of a type II OPO 
below threshold, paying special attention to the effects 
of walk-off. We have analyzed the effects of selecting dif- 
ferent polarization components: when walk-off vanishes 
or in the far field region not affected by walk-off (k y = 0), 
there is an almost complete suppression of noise in the 
proper quadratures difference of any orthogonal polariza- 
tion components of the critical modes (Sect. Ill A). Se- 
lecting non-orthogonal polarization components the cor- 
relations are reduced, vanishing for parallel polarizations. 
Walk-off strongly influences the strength of correlations. 
First, the variance of the quadratures difference of or- 
thogonal polarization components depends on the refer- 
ence polarization: the best EPR entanglement conditions 
are fulfilled when the most intense polarization compo- 
nents are locally selected ("vertical bright" scheme). If 
the selected polarizations are not the most intense ( "ver- 
tical dark" scheme) we still find EPR entanglement, but 
correlations are reduced and there is also a rotation of the 
quadrature angle giving the best squeezing (Sect.IIIB). 

Our study of EPR quadrature correlations identifies 
how these correlations depend on the polarization state. 
We have further investigated non classical polarization 
properties in terms of the Stokes operators. The proper- 
ties of the Stokes parameters in a single point of the far 
field do not show any non classical behavior (Sect. IV). 
For k v = 0, where there are no walk-off effects, we have 
shown that the average of the Stokes vector vanishes, 
and all the Stokes operators can be measured with per- 
fect precision. All Stokes operators are very noisy (above 
the level of coherent states) and the fluctuations are not 
sensitive to polarization optical elements: in fact the field 
is completely unpolarized, i.e. there is no 'hidden' polar- 
ization. Quantum effects are observed when considering 
polarization correlations between two symmetric points 
of the far field. Still in the direction orthogonal to the 
walk-off (k y = 0) we show perfect entanglement of all 
the Stokes operators measured in symmetric FF regions. 
This result is independent on the distance to the thresh- 
old. These results for k y = would apply in all the FF 
for vanishing walk-off. When there is walk-off and for 
k y 7^ the entanglement in the second and third Stokes 
operators is lost, but for So and Si there is still perfect 
correlation between two symmetric points of the FF, re- 
flecting the twin photon process emission. 
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FIG. 1. Far field intensity Eq. (19) for A = 0.99. Space 
is scaled with diffraction strenght ai = a,2 and time is scaled 
with cavity decay 71 = 72 [32], as reported in Ref. [13,15]. 
Other parameters are Ai = A2 = —0.25, pi — 0, p2 = 1. The 
same parameters will be used for all Figures, unless another 
choice is specified. The lower (upper) ring corresponds to the 
intensity of the field 1 (2). 
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FIG. 2. FF signal (black) and idler (grey) rings and de- 
tection scheme not influenced by the walk-off: The x and y 
polarizations are detected at the points ztfcn where the rings 
intersect (square symbols). 
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FIG. 3. FF signal and idler rings and detection scheme in- 
fluenced by the walk-off: a) The x and y polarized fields are 
detected in the points ±kv indicated by the square symbols, 
on the intense x polarized (black) and y polarized (grey) rings 
(bright detection), b) In the same points ±fcy (square sym- 
bols) the orthogonal polarizations y and x are detected (dark 
detection). 
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FIG. 5. V g (±k H ,Lu; $)/cr for * + *' = and g = 1 (dark 
continuous line), for VE' + = and g — g (light continuous 
line), and for ^ + <3/' ~ 0.06rad and g — 1 (dashed line). 



+ 



1.50 
0.75 

0.00 

-0.75 

1.50 
-3.24 




0. 1 



0.00 

CO 



3.24 



FIG. 4. V g =i{±k H ,u;$)/o for $ = (tt, 0, *, 0, tt/2, *') as 
a function of ^ + \P' (rad) and the frequency u). EPR entan- 
glement is obtained for values less than 1 (dark line). Values 
above 1 outside this line are not displayed. 
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FIG. 6. Same as in Fig. 4, but for the choice g = 
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FIG. 7. V fl= i(±fe v -,w;<E>)/cr for $ = (tt, tt/2, *, 0, n, *') as 
a function of if+^I/' (rad) and the frequency u>. The detection 
scheme is the one shown in Fig. 3a. EPR entanglement is 
obtained for values less than 1 (dark line). Values above 1 
outside this line are not displayed. 



12 



1.5 




3.0 



FIG. 8. Dependence of V g (±fcy, uj; $)/<t on the frequency 
uj, for \]> + = and g — 1 (dark continuous line), 
for + = and g — g (light continuous line), for 
$ + ~ 0.06rad and <j = 1 (dark dashed line), and for 
* + ~ 0.06rad and g = g (light dashed line). 
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FIG. 9. V g=g -(±kv,uj;$)/a as in Fig. 7, but for the choice 
9 = 9- 




FIG. 10. V s= i(±/W,^;<E>)/cr for $ = (vr, 0, *, 0, tt/2, *')> 
as a function of + (rad) and frequency w. The detection 
scheme is the one shown in Fig. 3b. 
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FIG. 11. Dependence of V g {±kv, w; 5>)/<7 on the frequency 
uj, for $ + = $ and g = 1 (dark continuous line), for 
$ + = $ and g — g (light continuous line). ^ ~ 0.81rad. 
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FIG. 13. Spectral variance Vg as a function of frequency, 
for the detection schemes represented in Fig. 2 (light line), 
Fig. 3a (dark continuous line), and Fig. 3b (dark dotted line). 
The angles VP + <&' = 0, 0, 0.81 rad are selected respectively 
for each detection scheme, a) walk-off p2 = 1; b) walk-off 
P 2 = 1.5. 



FIG. 14. Stationary average of the Stokes operator (Si(k)) 
in the FF (Eq.(52)). 
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FIG. 15. Second order polarization degree P%(k,i) defined 
in Eq. (54). The two circles of double continuous line show 
the signal and idler maxima intensities. 
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FIG. 16. Normal ordered variances Fi(fc,fc,0) and 
r2,3(fc,fe,0) normalized to the shot noise N = (So(k, t))a. 
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FIG. 17. D2(k,—k,0) normal ordered normalized to the 
shot noise. The white contour line shows the boundaries be- 
tween classical and quantum statistics. The two circles of 
double continuous line show the signal and idler maxima in- 
tensities. 
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FIG. 18. The same as in Fig. 17, but with walk-off p 2 = 0.5. 
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